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Spontaneous formation of permanent shear bands in a mesoscopic 
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This study proposes a coherent scenario of the formation of permanent shear bands in the flow of yield stress materials. It is 
a well accepted point of view that flow in disordered media is occurring via local plastic events, corresponding to small size 
rearrangements, that yield a long range stress redistribution over the system. Within a minimalistic mesoscopic model that 
incorporates these local dynamics, we study the spatial organisation of the local plastic events. The most important parameter in 
this study is the typical restructuring time needed to regain the original structure after a local rearrangement. In agreement with 
a recent mean field study [Coussot et ai, Eur. Phys. J. E, 2010, 33, 183] we observe a spontaneous formation of permanent shear 
bands, when this restructuring time is large compared to the typical stress release time in a rearrangement. The bands consist 
of a large number of plastic events within a solid region that remains elastic. This heterogeneous flow behaviour is different in 
nature from the transient dynamical heterogeneities that one observes in the small shear rate limit in flow without shear-banding 
[Martens et ai, Phys. Rev. Lett. , 20 1 1 , 106, 1 5600 1 ] . We analyse in detail the dependence of the shear bands on system size, shear 
rate and restructuring time. Further we rationalise the scenario within a mean field version of the spatial model, that produces a 
non monotonous flow curve for large restructuring times. This explains the instability of the homogeneous flow below a critical 
shear rate, that corresponds to the minimum of the curve. Our study therefore strongly supports the idea that the characteristic 
time scales involved in the local dynamics are at the physical origin of permanent shear bands. 
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1 Introduction 

Materials that only start to flow beyond a critical yield stress 
are ubiquitous in nature and serve in a wide range of applica- 
tions in industry and every day life. Famous examples in soft 
matter are foams, colloids, dense emulsions, gels, and pastes, 
but also granular matter and hard amorphous materials like 
metallic or polymer glasses fall into this category. At the mi- 
croscopic scale all of these materials have of course very dif- 
ferent properties and can barely be described by one single 
model. However on a more coarse grained level, when for 
example considering rheological or plasticity features, these 
materials show very similar properties. 

The steady state rheology in experiments with imposed 
shear rate y yields for example similar flow curves for a va- 
riety of complex fluids, that can roughly be categorised in two 
classes. The first class of yield stress materials will produce 
rheological curves that can be well fitted by generalised Her- 
schel Bulkley laws, a = oy + ay". In this description a is 

the average steady state shear stress and Oy is the dynami- 
cal yield stress, corresponding to the low shear rate limit. The 
prefactor a and the exponent n are material dependent fitting 
parameters. Such a law describes a monotonous behaviour 
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of the mean value of the stress as a function of the imposed 
shear rate. The second class of materials is characterized by 
flow curves that exhibit a minimum, indicating an instability 
of homogeneous flow in the region where da/dj < 0, and 
leading to heterogeneous flow in the form of shear bands 1 . 
These two classes have been distinguished in several recent 
works as "simple" fluids without shear banding behaviour and 
"thixotropic" fluids, that phase separate in a solid and a fluid 
phased - —. 

It is generally agreed that, as soon as the flow becomes 
heterogeneous, a description of the rheological behaviour in 
terms of a flow curve is insufficient. This has been illustrated 
e.g. in the experimental work of Goyon et al 6,7 , who showed 
that a single flow curve was not able to account for the flow 
profile of a colloidal paste in a microfluidic channel. Such a 
size dependence of the flow profile hints to the existence of 
length scales that are not part of the simple rheological de- 
scription, and several models have recently been put forward 
to account for the existence of such material specific length 
scales, either phenomenologically 89 or by an approximate 
coarse graining of a mesoscopic description of plastic flow 10 . 
The general structure couples a local rheological model with 
a diffusion equation that governs the behaviour of a "fluidity" 
parameter. As a result, coexistence between flowing regions 
with high fluidity, separated by an interface from regions with 
low plastic activity, is possible. 

The importance, in terms of strain localisation, of some 
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kind of local softening mechanism, has been pointed out by 
a number of groups (see e.g. reference^ for a discussion). In 
mesoscopic models such as the one we will be discussing be- 
low (see section 0, softening can arise from an initial ageing 
process, which results in a distribution of local yield stress 
that differs from the one in a "rejuvenated" system under 
flow—. This type of mechanism, that leads to transient (al- 
beit long lived) localisation of strain, was also studied within 
a "soft glassy rheology" context in reference 13 . The creation 
of shear bands due to structural relaxation (aging) has also 
been adressed in a detailed mesoscopic study by Jagl a 14 ' 15 . 

A different mechanism, which has been studied at a mean 
field level by Coussot and Ovarlez 16 , involves a local restruc- 
turing time that can be large enough to produce effectively a 
significant local softening after the system has undergone a lo- 
cal plastic deformation. Similar ideas, introducing directly a 
local weakening of the yield stress that persists over the dura- 
tion of avalanches after an initial slip event has occurred, were 
studied at the mean field level by Dahmen and co-workers 17 
and by Mansard et al. on the basis of a coarse-grained elasto- 
plastic model 18 . It can also be argued that a large restructuring 
time is, qualitatively, similar to the existence of an efficient re- 
laxation mechanism of the type described in reference s 14 i 15 . 
In both cases, the local stress will have time to relax signif- 
icantly before the element is loaded again. The purpose of 
the study reported in this pa- per is to investigate, beyond the 
mean field level, the transient softening mechanism proposed 
by Coussot and Ovarlez^. 

In the following sections we will show that a long local re- 
structuring time indeed gives rise to the formation of perma- 
nent shear bands. We will study the nature of these bands and 
show that this heterogeneous flow is different in nature with 
respect to the formation of dynamical heterogeneities in the 
simple flow at low shear rates 19 ' 2 ^. In section |2] we introduce 
the model we consider to describe the local dynamics on a 
mesoscopic scale. The rheological curves produced by these 
dynamics are compared to a mean field version of the model 
in section [3] In this section we motivate the importance of 
the local restructuring time as model parameter. In section |4] 
we study the spatial organisation of the plasticity for large re- 
structuring times as a function of the spatial form of the stress 
propagator, as well as the evolution in time of the shear band- 
ing instability. We show that the plasticity organises sponta- 
neously in bands, that persist in the long time limit. The long 
time characteristics of the dynamics are detailed in section [5] 
We show a dependence of the width and the density of the 
bands on the system size, on the typical restructuring time and 
on the imposed shear rate. In section|6]we summarize the most 
important results of our study and give a tentative explanation 
for the physical origin of the shear band scenario. 



2 Model 

To investigate in detail the influence of the different time- 
scales in the microscopic dynamics we will study a very sim- 
ple spatial (as opposed to mean field) model, that accounts 
for the elasto-plastic behaviour of the medium. The model is 
essentially the one introduced by Picard and coworker s 21 ' 22 . 
The medium is described by a set of elasto-plastic elements 
that occupy the nodes of a regular square lattice. To model 
a yield stress material under steady shear with a fixed strain 
rate y, we start from a simple Maxwell-like description of the 
stress dynamics associated with a given element: 

d t 0{t)=pLj-2pL£V\t) (1) 

where d t denotes the partial time derivative, a(t) is the (non 
tensorial) stress, jx is the shear modulus, and £ pl (f) accounts 
for the change in the strain due to local yielding. To take into 
account the long range effects of the plastic events, first pre- 
dicted by Argon 23 and later evidenced in molecular dynam- 
ics simulations^ 4 - and experiments on colloids 25 we change 
the above homogeneous equation to a field description with a 
stress propagator G(r — r') accounting for the stress redistri- 
bution due to a plastic event, 

d,a{r,t) =fiy+2ii J dr'G(r -r')e p \r\t) (2) 

To describe the deformation due to the plastic events, we 
model a Maxwellian viscoelastic like relaxation of the mate- 
rial in the plastic state, i.e. 

e p W) - 2^B(r,*Mr,0 . (3) 

where T is the typical time for the stress release in the plastic 
phase and n(r,f) a local state variable. In the following we 
refer to n(r,f) as local activity. We define n(r,f ) = indicating 
the absence of a plastic event and n(r, t ) — 1 if the local region 
is in the plastic phase. 

To incorporate the idea of a yield stress, we define the fol- 
lowing stochastic dynamics for the state variable «(r,f) (see 
left panel of Fig. [TJ. If the local stress exceeds a threshold 
value cr(r,f) > a y there is a finite probability to yield locally. 
The yielding rate is given by 1 /x ? \. Once yielded (n(r,f ) = 1) 
stress is redistributed using an appropriate stress propagator 
and locally the stress is relaxing during a typical restructuring 
time T e i. The rate to re-establish the elastic state (n(r, t ) = 0) is 
given by 1 / x s \. Note that in former studies of this model both 
rates have been set equal to the elementary relaxation time 22 . 
In this work we focus on the dependence of the dynamics on 
the ratio of these time scales. 

Approximating the rearrangements corresponding to the lo- 
cal plastic events by a force quadrupole (Eshelby problem), 
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one expects a four-fold quadrupolar symmetry for the inho- 
mogeneous part of the stress propagator 26 (see right panel of 
Fig. [TJ. The form of this propagator for an infinite two dimen- 
sional medium reads 

G H {r,0) = -i-rcos(40) . (4) 




Fig. 1 Local stress dynamics: Left: Schematic representation of the 
local stress of an isolated site as a function of time. When the stress 
is larger than the threshold o" v the site yields with a probability given 
by the plastic rate While the site is active the stress drops 

exponentially during a typical restructuring time T e i until it becomes 
elastic again. The site is supposed to be isolated, to simplify the 
picture. In the completely connected model described above we 
would expect stress fluctuations due to the mechanical noise coming 
from the surrounding sites. Right: Stress change due to a single 
plastic event at the origin. The stress change at the origin is 
artificially set to zero to account for the visibility of the nearest 
neighbour changes. In reality the change at the origin is strongly 
negative. 

In the following we will restrict our interest to athermal sys- 
tems in a planar geometry. Further, to simplify the problem we 
consider an incompressible medium and assume that the mi- 
croscopic geometry of the plastic events is the same as the one 
of the macroscopic shear. This simplification permits to use 
the above introduced scalar description of the stress. Note that 
convection effects on the stress and the activity have been ne- 
glected to keep the model as simple as possible. The lack of 
convection leads basically to an additional x-y symmetry for 
the spatial arrangements of the plastic events. For the shear 
bands this means that there is no preferred direction of the 
flow, such that they will be oriented with equal probability in 
the x or y direction. 

To simplify, we rewrite Eq. (O in dimensionless quantities, 
where the time t is measured in units of the elementary relax- 
ation time T, the stress a in units of the local yield stress a y 
and the strain rate y in units of the critical value % = Oy/tpL. 
% indicates the change from a Newtonian to a non Newtonian 
flow behaviour. To concentrate on the bulk dynamics we im- 
plement periodic boundary conditions and we integrate Eq. [2] 
numerically following the approach of Picard et al. 21 ' 22 . We 
encode the local stress values Oy (i and j being the discretized 



space coordinates) on a square lattice of size N = L 2 , given in 
units of a 2 , the typical size of a plastic event. Further we cal- 
culate the convolution in the second term of Eq. [2] using fast 
Fourier transform methods. We approximate the finite size 
propagator corresponding to Eq. |4]in Fourier space as 

G( qi ,qj) = -4 - ' 1 (5) 
{Ti+TjY 

with wave vector q = (qi,qj) = ^r(i,j). The homogeneous 
part of the propagator is implemented as G(0,0) = — 1. The 
integration time step is taken to be small with respect to all 
other time scales, Af/r = 10~ 2 . 

The numerical implementation of the stochastic yielding 
dynamics is straightforward. To be precise, we summarize 
the above described dynamics for the activity field ny in the 
following expressions 

P ol (n i j(t + dt) = l\n i j(t) = 0;a l j> a y ) = — (6) 
P l o{n i j(t + dt)=0\n iJ (t) = l) = — , (7) 

^el 

where Poi(n(t + dt) = i\n(t) = 0;rj > oy) denotes the tran- 
sition probability of going from state n = to state n = 1 in 
a small time interval dt, given that the local stress value ex- 
ceeds the yield stress. Pio(n(t + dt) — 0|«(f) = 1) denotes the 
corresponding reverse transition probability. 

After the rescaling we obtain a minimal model describing 
the dynamics of two fields, namely the stress field ay and the 
activity field ny with only a few parameters, that all have a 
physical motivation. The system size N and the shear rate y 
are macroscopic controllable quantities. Further there are the 
different time-scales involved in the microscopic dynamics, 
T p i and T e i, describing the typical time to yield and the typical 
restructuring time respectively. The yielding time T p i accounts 
for the disorder in the system. Defining a probability to yield 
above the globally defined yield stress is equivalent to intro- 
duce random barriers in the system, that change once a site 
has yielded. This is similar to the definition of the soft glassy 
rheology model 27 , but with a strain rate dependent probability 
distribution for the barrier heights. The simulations presented 
in this paper have been performed for a value of the typical 
time to yield equal to the elementary relaxation time (T p i — 1). 

The most interesting control parameter in this study is how- 
ever the restructuring time T e i. This parameter determines how 
long the system is locally in a plastic state after a yielding 
event. The larger T e i the longer we relax locally the stress dur- 
ing a plastic event and redistribute it over the system. In a 
more physical interpretation we expect this parameter to cor- 
respond to the time the material needs to regain its original 
structure, after links are broken trough the local rearrange- 
ments in the plastic event. This interpretation is similar to 
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what has been recently proposed by Coussot et al. 16 within a 
mean field model for shear banding behaviour. In the follow- 
ing sections we will show that this parameter plays an essential 
role in the formation of a shear band structure. 

3 Rheology 

Before going into the study of the spatial organisation of the 
plasticity within our model, we shall study some macroscopic 
features, such as the flow rheology. To do so we derived an 
analytically solvable mean field description for our model, in 
which we describe the effective dynamics on a single site. We 
rewrite the stress dynamics of the full model as 



dtOijit) = tff (t) + G(0,0)n i7 (0M ? ) 



(8) 



where G(0,0) < is the value of the stress propagator at the 
origin; <7; ; -(f) denotes the shear stress at point {i,j} and 



trim 



(9) 



If we assume yfj 1 to be homogeneous in space and if we 
neglect any possible fluctuations around its mean value f we 
obtain the following effective dynamics for the stress evolu- 
tion on one site 



d t a{t) = t-gn(t)a(t) 



(10) 



with g = — G(0,0) > the absolute value of the stress propa- 
gator at the origin. This value is system size dependent in the 
fully spatial model, but approximately constant in the limit of 
large systems. In the range studied within this paper g is large 
with respect to the shear rates of interest {g rs 0.57). The ac- 
tivity n(f) obeys the same local stochastic rules as in the full 
spatial dynamics, given in equation (|6). The advantage of this 
mean field model is that we can extract analytically the expres- 
sion for the flow curve, that is the average steady state stress 
as a function of the effective shear rate f . 

To derive an approximate solution for the mean field equa- 
tion ( [Tol l, we assume f < g and T e i sufficiently large to decor- 
relate the typical stress values for the change of activity n. We 
have the following two equations that describe the average dy- 
namics in the inactive (n = 0) and in the active (n = 1) part of 
the dynamics: 



« = 0: (7(f) = G- 



Yt 



n=l: a(f) = - +A(f, Tpi)e 



-gt 



(ID 
(12) 



where the bar indicates an ensemble average. Further (7 is 
given by the typical stress value at which the site becomes 
elastic 



a 



8 1+gTel 



(13) 



and the prefactor A(f , T p i) = (7+ — t/g is fixed by 

CT+ = l+fT p i, (14) 

the typical stress value at which the site yields. 

To calculate the time averaged stress (o)t, we now aim for 
the analytical expressions for the typical times spent in the 
two different activity states and the associated average stress 
values 



(<*>» = 



Tin ( (7in ) r + T ac t ( "act ) 1 



(15) 



by 



Tin T" T ac t 

The time average stress value in the inactive phase is given 
(a m ) t = — I da+P+{a+) j da-P-(a-) x 



(7+ + G \ / (7+ - (7 



(16) 



where P + and P denote the probability distributions of (7_ 
and o+ respectively. The typical corresponding time, spent in 
the inactive phase simply yields T; n = ( (7+ — (7 ) /f = T p i + 
(1 — (7_)/r. Therefore we obtain for the average stress in the 
inactive state 



with 



2((7+-(7_) 



l+2fT pl a 



rv 2 _fA(f,T pl ) 



(17) 
(18) 

, - , d9) 

8/ 8 1+gTel 

^-2(r/g)<r + + (f/g) 2 

l+2gT d 

The last missing information is the average stress value in 
the active state, (7 act which can be calculated by 



(<7act)r = -f 



' dxe z l^ f dt'(- 
JO \g 



f | A(f,T pl ) 

8 1+gTel 



(20) 



This means that the average stress in the active state is equal 
to the typical stress value to become elastic ((7 act ) r = ~o . The 
typical time spent in the active state is simply given by T act = 

Tel. 

Using the analytical expressions in the former paragraph we 
can express the steady state average stress as a function of 
strain rate. In Figure [2] we show the steady state solution of 
these mean field equations for various values of the restructur- 
ing time Tei. We find that above a critical value of this restruc- 
turing time Tgj the mean field solution exhibits a minimum. 
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This means for a given value of the average stress (a) we find 
more than one value of the effective shear rate, which may give 
rise to phase coexistence in the spatial model. From the above 
equations the critical value for the restructuring time can be 
derived through limp^ d (a) /dT = 0, yielding 



C, = 



2g 



(21) 



For a typical yielding rate = 1 and an absolute value of 
the propagator at the origin of g = 0.57 this yields a critical 
restructuring time tS « 3.3. For restructuring times above this 
critical value we obtain a minimum in the flow curve within 
the presented mean field description. There is a critical strain 
rate F c below which the rheological curve exhibits a negative 
slope. The value of this critical shear rate will both depend 
on the restructuring time T e i and on the yielding time T p i (see 
right panel of Figure |2). 
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Fig. 2 Rheological curves: Left: Shown are the numerically 
measured values of the steady state stress (<j) of the mean field 
model as a function of strain rate f for different values of the 
restructuring time (T e ] = 1,2.5,5,7.5, 10) for T p i = 1. This is 
compared to the analytically derived expressions for the stress 
dependence on the strain rate for the same restructuring times (solid 
lines). Right: Strain rate T c corresponding to the minima of the 
mean field flow curves in the left panel as a function of restructuring 
time T e ] for different values of the yielding rate 
(T pl = 1.0,2.0,3.0,4.0,5.0). 

Another quantity of interest, that can be easily derived in 
this mean field approximation, is the average steady state ac- 
tivity (n)t as a function of effective shear rate T, given by 
(n) t = T act /(T act + Ti n ). This yields 



(n)t 



g(?el + T p l + l/r) 



(22) 



An alternative way to derive the average activity is to consider 
the steady state state expression (d,a) r = and thus (no) t = 
T/g, which can be rewritten as (n)t<J a a =r/g, recovering the 
same expression I 



The main difference between our mean field model and the 
former introduced model of Coussot et al. is the existence of 
a third time scale, which is the typical rate to yield T„/. In the 
model of Coussot et al, the flow behaviour turns from that of 
a simple yield stress fluid to a shear-banding material when 
the ratio of a characteristic relaxation time of the system to a 
restructuring time becomes smaller than one 16 . Note that in 
the limit of T p i = and a value of g = 1, we recapture this 
result within our mean field model (see Eq. 1211 1. The flow 
curve however is different from the simple mean field result in 
the work of Coussot and coworkers, due to the more complex 
relaxation dynamics (exponential decay of the local stress). 

4 Band formation 

From the above observed minimum in the rheological curve 
we expect that the activity localises in some specific regions. 
But are the plastic events really concentrated in shear bands? 
One way to get an idea about the spatial organisation of the 
plastic events is to calculate numerically the accumulated ac- 
tivity map, that is the local number of plastic events in a given 
strain window. It turns out that the specific form of the stress 
propagator plays an important role for the spatial distribution 
of the plastic events. When the original long range propaga- 
tor (given in Eq|4|i is used to redistribute the stresses after a 
plastic event occurred, we indeed obtain a shear band struc- 
ture (see right panel of Figf3j. However if we use a simplified 
short range propagator with the same homogeneous part but 
an inhomogeneous stress redistribution, given by the follow- 
ing table: 



-1/4 


1/2 


-1/4 


1/2 


-1 


1/2 


-1/4 


1/2 


-1/4 



we obtain active regions of typical size, but of a shape that 
is very different from the band geometry (see right panel of 
Fig- ID. This means that the long range character of the prop- 
agator enhances the shear banding effect. This can be un- 
derstood from the fact that, with a short range propagator, a 
plastically active region of arbitrary shape can coexist with an 
inactive region, with a smooth interface separating the two re- 
gions. For the long range elastic propagator, on the other hand, 
a uniform line of slip events is the only type of plastically ac- 
tive region that does not give rise to any deformation in the 
surrounding medium. 

To test whether these bands are only a transient phenom- 
ena leading to a homogeneous state in the long time limit, we 
performed the following simulation. At the beginning we ap- 
ply a high strain rate 7=1 where the dynamics is homoge- 
neous, that is the probability of producing a plastic event in 
the stationary state is uniform in space. Then we perform a 
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Fig. 3 Cumulated activity: Local number of plastic events in a 
given strain interval Ay = 100 in a system of size N = 2 16 with a 
restructuring time T e j = 10.0. Left: stress is redistributed using the 
standard long range stress propagator, given in Eq. I0. Right: stress 
dynamics are governed by an alternative short range stress 
propagator (see text). 



rapid quench of the shear rate into the unstable regime with 
y = 0.01. This allows to start from a homogeneous condi- 
tion. Moreover due to the periodic boundary conditions and 
homogeneous applied shear rate, there is no way to break the 
homogeneity except for the appearance of instabilities due to 
collective effects in the stress dynamics. And indeed we ob- 
served the appearance of thin and short shear bands and an 
ensuing coarsening that lead to longer and wider bands, even- 
tually reaching the system size (see Fig. @). In the long time 
limit one system size spanning band is formed, that has a char- 
acteristic width and density. This band can drift in perpendic- 
ular direction of its orientation. 

Note that for very small system sizes, strong fluctuations 
can lead to a renucleation of the band in another region. But 
even in this regime the structure factor remains constant and 
we obtain a permanent instability. For very large system sizes 
we observe a transient regime in which several bands can co- 
exist during a finite time. However, for the parameter range 
studied here we always obtain a steady state (with a time inde- 
pendent static structure factor) producing one single band that 
only moves perpendicular to its orientation. The system size 
is chosen sufficiently large to ensure that fluctuations will not 
lead to a renuclation of the band elsewhere in the system. 

5 Long time characteristics of the shear bands 

In this section we will discuss the long time characteristics 
of the formed bands. As we showed in the former section 
the spatial organisation of the activity in the long time limit 
corresponds to one single band of a given width, density and 
interface width. These characteristics of the band will depend 
on the model parameters, that is system size N, restructuring 
time T e i an d strain rate y. 

But before studying in detail these dependences we need to 
give a definition of what we want to call a permanent shear 
band. To obtain a characteristic band profile we need to aver- 
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f Ay 0.5 


Ay = 1 .0 
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Fig. 4 Quench dynamics: Instantaneous activity maps nu, with 
i, j € [0, L] at different strain intervals Ay after a rapid quench of the 
shear rate from y = 1.0 to y = 10~ 2 for a system size N = 2 16 and a 
restructuring time of T e ] = 10. Active regions are indicated in red 
(ny = 1), passive regions in white (n/j = 0). 



age over a certain strain window Ay w . This strain window must 
be large enough to capture the dynamically built shear band, 
but small with respect to the typical strain that causes the band 
to change its position. In the left panel of Fig. [5] we show the 
band profile as a function of different strain windows, used for 
the averaging. Indeed, we find that it is possible to obtain a 
band profile that is well defined for a range of averaging strain 
windows, as long as these are small with respect to the strain 
window that leads to a significant band displacement (perpen- 
dicular to the band orientation). 

To obtain the long time characteristics of a band we use a 
gliding window average in the steady state, where we trans- 
late the band, obtained in each window, to the same position 
(centre of the system) for the averaging. In the right panel of 
Fig. [5] we show such a band profile for given model parame- 
ters that yield inhomogeneous flow. Indeed we find a region 
with a strong activity that we call shear band and regions with 
almost no plastic events, that we will call solid phase in the 
following. 

But how does the profile of the average activity (n) compare 
to the profile of the average effective shear rate? To obtain this 
relation we study the steady state dynamics where d t (a) = 0, 
which yields, using relation i 

1 



Y+jjCE, Gii-i'J-fni'/Oitj^t = 
u'jf 



(no) 



y. (23) 



This reflects the fact that the space averaged change of plastic 
strain equals the globally applied shear rate, leading to a non- 
equilibrium steady state. 

Further it turns out that for homogeneous flow, T e i < tS ac- 
tivity and stress are uncorrelated, such that (n)(a) = j (see 
left panel of Figure [6] This is in agreement with the KEP 
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Fig. 5 Band characteristics Left: Average activity (ft) as a function 
of space (coordinate perpendicular to the band orientation). The 
average is taken over the second space coordinate and over a given 
strain window Ay,, (10, 100, 200, 300, 400, 500, 600, 700, 800, 
900). The inset shows the rescaled position of the band x p /L as a 
function of strain from band profiles obtained through averaging 
over a strain window of Ay H , = 10. Right: Average local effective 
shear rate {y e ff)and rescaled average local activity (n) jC as a 
function of space for restructuring time T e ] = 10, shear rate 
y= 10~ 16 and system size N = 2 16 . The numerically measured 
value of the rescaling factor is C = 0. 155 (see text). 

model assumption—, if we interpret the local activity as flu- 
idity. However in the inhomogeneous flow regime there are 
strong negative correlations between activity and stress (see 
right panel of Figure [6]). This is intuitive since an active site 
that remains active for a long time will lead to a strong local 
stress relaxation. 

To get an estimate of the relation between the effective shear 
rate and the activity, we study equation ([H) in the steady state 
regime where we find 

(ff)t = g(mjOij)t 

~ g( n ij)t Oact , (24) 

where a act indicates the time averaged stress value on an active 
site. We find that the local average effective shear rate is linear 
in the local average activity, (yfj)t ~ C(«, / ) r , with C — g(7 act . 
As mentioned earlier the value of g depends on the systems 
size. For example for the parameters N = 2 16 , y= 10~ 16 
and T e i = 10 the proportionality constant yields C sa 0.155, 
see Fig|5] 

It immediately follows that in the regions where (tin) ~ 
we also obtain a vanishing small effective local shear rate. We 
will refer to these regions as solid regions. On the contrary 
in the regime with high activity we obtain a finite effective 
shear rate corresponding to a shear band (see right panel of 
Fig 0. In the following we will use the activity rather than 
the measurement of the effective shear rate for the definition 
of the shear bands. 

Let us now study in detail the long time characteristics of 
the band as a function of system size, shear rate and restruc- 




Fig. 6 Flow and activity curves: Left: Shown are the average stress 
(cr), the average activity (ft) and the product of both (n)(a) as a 
function of strain rate y, for a system size N = 2 16 and a value of the 
restructuring time T e \ = 1 (homogeneous flow). Right: Shown are 
the average stress (o"), the average activity (n), the product of both 
(;i)(o") and the average stress weighted by the activity of the 
corresponding site (nd) as a function of strain rate y, for a system 
size N = 2 16 and a value of the restructuring time T e ] = 10 
(inhomogeneous flow). 

turing time. In the left panel of Fig. [7] we display this band 
profile as a function of system size for T e i = 10 and f= 10~ 2 . 
It turns out that the width of the band is linear in system size 
and the width of its interface £, has no significant dependence 
of the size of the system (see right panel of Fig. [7). 




Fig. 7 System size dependence of the band profile: Left: System 
size dependence of the band profile. Shown is the average activity 
(ft) in the steady state as a function of space (coordinate always 
perpendicular to the band orientation) rescaled with system length L 
for different values of the system size N = L 2 . The averaging strain 
window is Ay tt . = 10 and we average over a total strain of Ay = 1000. 
Right: Space-derivative of the average activity (n) as a function of 
space for different values of the system size. The position of the 
band is shifted, such that the beginning of the band is in the same 
position. The length t, is a measure of the width of the interface. 

As for what is the dependence of the band characteristics on 
the shear rate, we find that the maximum (or peak) activity of 
the band does not depend on y. The peak activity h is here de- 
fined as the difference between the minimal and the maximal 
value of the average activity. The average effective shear rate 
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in the shear band takes the value where the average stress 
in the corresponding effective flow curve is minimal (see Fig- 
ure|2]l. This fixes the average stress in the solid phase as well. 
It is only the relative width of the band w that changes with 
y. For large y the dependence is linear. This allows for an 
extrapolation, yielding the critical shear rate % w 5.5 x 10~ 2 
at which the width of the band reaches the system size. This 
is the point of the crossover from a heterogeneous flow to a 
homogeneous flow without any shear bands (see left panels of 
Figure [8] 




Fig. 8 The flow transition: Upper left: Band profile of the single 
band forming in the long time limit as a function of the shear rate y 
(1CT 2 ' 2 , 1(T 2 - , 1CT L8 , 1(T L6 , 1(T L4 ) for a restructuring time 
T e i = 10. Upper right: Band profile as a function of the 
restructuring time T e , (5.5, 6.5, 7.5, 10, 12.5, 15, 17.5, 20) for a 
shear rate y = 10~ 2 . Lower left: Width w (rescaled with system 
length) as a function of y for T e ] =10. The dashed line is a linear fit 
of the large shear rate region and the black circle indicates the point 
where we expect the transition to homogeneous flow (the width of 
the band reaches the system size at y c a> 5.5 x 10~ 2 ). The dotted line 
separates the region of homogeneous from the region of 
inhomogeneous flow. Lower right: Inverse peak activity hT x of the 
band profile as a function of the restructuring time T e ] for y = 10~ 2 . 
The dashed line is a fit (see text) and the dotted line indicates the 
fitted critical restructuring time for the transition from homogeneous 
to inhomogeneous flow behaviour at xS ~ 4.3 (see text). 

When changing the value of the restructuring time, we ob- 
serve a different phenomenon. For large T e i the width of the 
band is fixed, but the peak activity h of the band grows with 
increasing T e i. On the other hand below a critical value of the 



restructuring time the band vanishes and the flow becomes 
homogeneous. To estimate this crossover time, we use that 
in the steady state (n)a act equals y. In the parameter range 
where the width is fixed the peak activity must compensate 
the decreasing value for the average stress on the active sites 
(7 act due to the long time of local relaxation. A simple fitting 
function for the inverse peak activity is 

h-' = ^-+h-J (25) 

Tel ~ ^ 

with the fitting parameters /i~ 1 > 1 , being the peak activity at 
large restructuring times, a{y), a shear rate dependent prefac- 
tor and T^, the critical restructuring time for which the system 
becomes homogeneous. With this fitting procedure we find a 
value of the critical restructuring time of sa 4.3, which is in 
fair agreement with the findings within the mean field model, 
see Eq. (l2Tb . Note, that for very large T e ] the time to form 
one single band grows strongly and becomes inaccessible. In- 
stead of one large band we observe several smaller coexisting 
bands. Thus, the above relations are only valid in a finite range 
of values for T e i. 

6 Discussion and conclusions 

We performed the first detailed spatial analysis of the forma- 
tion of shear bands within a minimalistic model, that bridges 
the microscopic dynamics to the macroscopic flow behaviour. 
The model for the yield stress flow invokes only few physi- 
cally relevant parameters from the microscopic and from the 
macroscopic scale. 

We find that a microscopic time-scale, that is linked to the 
restructuring time of the material, controls the crossover from 
homogeneous to heterogeneous flow, similar to what has re- 
cently been predicted within a mean field study—. If the time 
needed after a local plastic event, to regain the original struc- 
ture is large in comparison with the typical duration of local 
stress release, the system undergoes a phase separation below 
a critical shear rate. The separation gives rise to a strongly 
sheared phase of band geometry within a solid non flowing 
environment. We show that this band geometry is closely re- 
lated to long range character of the stress redistribution due to 
a local plastic event. 

After a quench from a large to a small shear rate we observe 
coarsening behaviour for the activity in the system, that yields 
in the long time limit a single strongly active band. The band 
characteristics are well defined and depend on the system size, 
the value of the shear rate and the typical local restructuring 
time. We study these dependences in detail and show how to 
get a prediction for the relevant critical values at which the 
crossover from simple to shear band flow occurs. 

These findings are compared to a mean-field version of the 
model that incorporates the influence of all other surrounding 
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sites into an effective shear rate. There are two important ques- 
tions to answer with respect to the spatial model to be able to 
make predictions about the phase coexistence. Let us assume 
that our dynamics produce two different phases in form of a 
shear band and inactive solid regions. Then we have to show 
that the average stress will be the same in the two phases. And 
indeed the stress dynamics are defined such that if we start 
from an homogeneous initial condition for the stress field, the 
dynamics will assure that the average stress in a line or in a 
column is equal to the global average stress. This is due to the 
specific form of the stress propagator, that is defined such that 



£g(u) 



1 

L 



(26) 



The change of the average stress on one line (horizontal or 
vertical) during one time step At turns out to be independent 
of the position of the line 



Y + \ L G (' ~~ i'lJ ~ i') n i'j' °i' fj At 



(Aoij}i = {f-(nijOij)ij)At 



(27) 



which is independent of j. The same argument is valid for 
summation over space coordinate i. If we start from a ho- 
mogeneous condition for the stress values (in the simulations 
Oy = 0, V/j) the average stress on a line will always be equal 
to the global average of the stress in the system. Thus, by 
construction the average stress in a shear band will equal the 
average stress in the inactive phase. The average stress can 
therefore be interpreted as intensive thermodynamic parame- 
ter that governs the phase coexistence. 

Further we must check whether the mean field description is 
relevant for the spatial model. More precisely we have to test 
whether fluctuations of the effective strain rate are really irrel- 
evant in the dynamics. We calculated numerically the steady 
state average stress as a function of the steady state average 
effective shear rate for a choice of the restructuring time for 
which we would expect the appearance of phase separation 
(see right panel of Figure [9). For large enough strain rates it 
turns out that the mean field model reproduces very well the 
dynamics of the spatial model. However for smaller strain 
rates there is a discrepancy between the two curves. This dif- 
ference results from the stronger impact of the fluctuations of 
the effective shear rate on the dynamics, yielding a shift of the 
minimum of the flow curve towards smaller effective strain 
rates. This means that the mean field model underestimates 
the parameter range of homogeneous flow. But qualitatively 
we observe very similar dynamics as predicted by the mean 
field model. For effective strain rates below a critical value 
%s(%) the flow becomes heterogeneous with two effective 
strain rates. This accounts for the existence of a frozen phase 
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Fig. 9 Phase separation scenario: Left: Sketch of a typical mean 
field flow curve that leads to a phase separation. The green line 
indicates all possible values for the stress and the blue line all 
possible values for the high shear rate phase. In red are indicated the 
dynamically chosen values for the phase coexistence, j e S = and 
Yeff{%) at equal stress {a) c , that correspond to the minimum of the 
flow curve. Right: Average steady state stress value {o) ss as a 
function of the effective strain rate %g for a value of the 
restructuring time T e i = 10 in the complete model with A' = 2 16 (*) 
in comparison with the flow curve of the mean field version of the 
model (x). The vertical dotted line indicates the critical value of the 
effective shear rate that separates the region with 

homogeneous flow from the region with shear band flow. The two 
values of y e ff in the heterogeneous regime have been measured 
separately in the two phases. The inset shows the width w of the 
band as a function of the strain window Ay, used for the averaging. 



with feS ~ and a strongly fluidized phase with y e ff « y e ff (%). 
As shown in the inset of the right panel of Figure [9] the ef- 
fective strain rate is linear in y for y > % and the effective 
strain rate in the shear band region is independent of y for 
y < %. Note that this scenario implies that the macroscopic 
flow curve, (a) as a function of y, is a monotonic curve, not 
easily to be distinguished from a flow curve that corresponds 
to homogenous flow. Only the spatial analysis and the corre- 
sponding mean field study reveal the inhomogeneous nature 
of the flow. 

In conclusion we can say that the mean field model pro- 
duces qualitatively a good description of dynamics of the fully 
spatial model (see left panel of Fig. [9] The qualitative be- 
haviour of the two flow curves are the same. Both exhibit a 
minimum in the average stress as a function of shear rate for 
large restructuring times. This minimum explains the phase 
separation since for a given average stress, that is homoge- 
neous in the system, two phases can coexist, one with zero 
effective shear rate and one with a finite effective shear rate. 
Hence, we provide a coherent scenario for the flow instability. 
Our results support strongly the idea that time-scales in the 
microscopic dynamics play an important role for the macro- 
scopic flow behaviour. 

Linking mesoscopic models such as the one studied above 
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with 'real world' systems is a difficult issue. One may proceed 
as Coussot and Ovarlez by infering from the consequences 
(namely the observed stability of bands) that some systems, 
such as attractive colloids, correspond to large T e i, while other, 
such as foams, correspond to small T e i. It is also interesting to 
speculate on the type of microscopic features that may give 
rise to a large restructuring time. Very generally, one can 
imagine that the restructuring time will be large in all cases 
where the shear can bring the system into a locally trapped 
configuration. The mechanical properties resulting from this 
metastable arrangement are thus expected to be significantly 
weaker than those of a relaxed configuration. Such a situation 
will be encountered for example in the case of strongly direc- 
tional interactions (e.g. patchy colloids or covalently bonded 
systems 28 ). Another potentially interesting situation is the 
case where the interaction is mediated by a different compo- 
nent which may have a long relaxation time, as e.g. is the 
case for capillary bridges in some emulsions or in wet gran- 
ular matter. One may also consider pair interactions between 
particles that involve an oscillatory potential and therefore a 
finite barrier between a 'bonded and a 'non bonded' state be- 
tween neighboring particles. Interestingly, such a potential has 
been used for a long time to describe the formation of a single 
component quasicrystal 29 , which is also a good glass former. 
It will be of interest to study the statistical properties of sim- 
ple flows in such models in order to obtain a mapping onto our 
more coarse grained description. 
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